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O ■ ABSTRACT 
(N 

, Aims. The goal of this work is to constrain the strength and structure of the magnetic field associated with the environment of the 

^ ■ radio source 3C 449, using observations of Faraday rotation, which we model with a structure function technique and by comparison 

' with numerical simulations. We assume that the magnetic field is a Gaussian, isotropic random variable and that it is embedded in the 

hot intra-group plasma surrounding the radio source. 

Methods. For this purpose, we present detailed rotation measure images for the polarized radio source 3C 449, previously observed 
with the Very Large Array at seven frequencies between 1.365 and 8.385 GHz. All of the observations are consistent with pure fore- 
ground Faraday rotation. We quantify the statistics of the magnetic-field fluctuations by deriving rotation measure structure functions, 
which we fit using models derived from theoretical power spectra. We quantify the errors due to sampling by making multiple two- 
I dimensional realizations of the best-fitting power spectrum. We also use depolarization measurements to estimate the minimum scale 

• . of the field variations. We then make three-dimensional models with a gas density distribution derived from X-ray observations and a 

i-^ ' random magnetic field with this power spectrum. By comparing our simulations with the observed Faraday rotation images, we can 

determine the strength of the magnetic field and its dependence on density, as well as the outer scale of magnetic turbulence. 
Results. Both rotation measure and depolarization data are consistent with a broken power-law magnetic-field power spectrum, with 
?^ , a break at about 1 1 kpc and slopes of 2.98 and 2.07 at smaller and larger scales respectively. The maximum and minimum scales of 

the fluctuations are ^65 and a;0.2kpc, respectively. The average magnetic field strength at the cluster centre is 3.5+1.2/jG, decreasing 
^ ' linearly with the gas density within a; 16 kpc of the nucleus. At larger distances, the dependence of field on density appears to flatten, 

but this may be an effect of errors in the density model. The magnetic field is not energetically important. 

Key words. - Magnetic fields - Polarization - galaxies: ISM-radio continuum: galaxies-X-rays: galaxies 

^ 1. Introduction Pollack et al. 2005; Govoni et al. 2006; Guidetti et al. 2008). 

QP , • ^ , , • , , , ■ , ■ , , The central magnetic field strengths deduced from these data are 

P Magnetic fields m the hot plasma associated with groups and ^^^^jj ^ ^ ^^^^^^ q .^^^ 

^ clusters of galaxies are poorly understood, but are thought to ^.^j^^^^ cool-core clusters (e.g. Taylor et al. 2002). The RM dis- 

O ' . • " '^^"l conduction (e.g. Balbus jributions of radio galaxies in both interacting and relaxed clus- 

O : 2000; Bogdanovic et al 2009) and influencing the dynamics ^^^^ generally patchy, indicating that cluster magnetic fields 

^ , of cavities formed by radio jets (e.g. Dursi & Pfrommer 2008; ^^^^ ^^^^^^j^^^ ^^^j^^ ^ Iq,^ 
O'Neill et al. 2009). The existence of magnetic fields can be 

demonstrated in several different ways (e.g. Carilli & Taylor Several studies of Abell clusters (Murgiaet al. 2004; Govoni 
V 2002; Govoni & Feretti 2004 and references therein). One of et al. 2006; Guidetti et al. 2008) have shown that detailed RM 
the ways of studying these fields is via the Faraday effect: ro- images of radio galaxies can be used to infer not only the 
a ' tation of the plane of linearly polarized radiation by a magne- sti-ength of the cluster magnetic field, but also its power spec- 
■ " " ' tized plasma. Synchrotron emission from radio sources (either trum. The analysis of Vogt & EnBlin (2003, 2005) suggests that 
behind or embedded within the group/cluster medium) can be the power spectrum has a power law form with the slope appro- 
used to probe the distribution of foreground Faraday rotation, priate for Kolmogorov turbulence and that the auto-correlation 
These can be combined with X-ray observations (which provide length of the magnetic field fluctuations is a few kpc. The de- 
the thermal gas density profile) to infer the strength and fluctua- duction of a Kolmogorov slope could be premature, however: 
tion properties of the magnetic field. there is a degeneracy between the slope and the outer scale 
Faraday rotation studies of clusters have been carried out us- which is difficult to resolve with current Faraday-rotation data 
ing both statistical samples of background radio sources (e.g. (Murgia et al. 2004; Guidetti et al. 2008; Laing et al. 2008). 
Lawler & Dennison 1982, Vallee et al. 1986, Kim et al. 1990, Indeed, Murgia et al. (2004) pointed out that shallower magnetic 
Kim et al. 1991, Clarke et al. 2001) or individual radio sources field power spectra are possible if the magnetic field fluctua- 
within the clusters (e.g. Taylor & Perley 1993; Feretti et al. 1995; tions have sti-ucture on scales of several tens of kpc. Recendy, 
Feretti et al. 1999a,b; Govoni et al. 2001; Eilek & Owen 2002; Guidetti et al. (2008) showed that a power-law power spectrum 

with a Kolmogorov slope and an abrupt long-wavelength cut- 

Send ojfprint requests to: dguidett@eso.org off at 35 kpc gave a very good fit to their Faraday rotation and 
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depolarization data for the radio galaxies in A2382, although a 
shallower slope extending to longer wavelengths was not ruled 
out. 

While most work until recently has been devoted to rich clus- 
ters of galaxies, little attention has been given in the literature 
to sparser environments, although similar physical processes are 
likely to be at work. Faraday-rotation fluctuations have previ- 
ously been detected in galaxy groups (e.g. Perley et al., 1984; 
Feretti et al. 1999a), but without deriving in detail the geometry 
and structure of the magnetic field. The first detailed work on 
galaxy groups was done by Laing et al. (2008), who analysed 
the radio emission of 3C 3 1 . 

They found that the three-dimensional magnetic-field power 
spectrum vv(/),defined in Sect. 16.11 might be described in terms 
of spatial frequency / by a broken power law w(f) oc Do/ 
with q - 1 1 /4j for / >0.062 arcsec (corresponding to a spa- 
tial scale of about 17kpc) and q=2.32 at lower frequencies, al- 
though a power spectrum with a slope of 2.39 and an abrupt 
cut-ofF at high frequencies could not be ruled out. Their results 
are qualitatively similar to those for sources in Abell clusters. 

Magnetic fields associated with galaxy groups deserve to 
be investigated in more detail, since their environments are 
more representative than those of rich clusters. Moreover, 
observations, analytical models and MHD simulations of galaxy 
clusters all suggest that the magnetic-field intensity should 
scale with the thermal gas density (e.g. Brunetti et al. 2001, 
Dolag 2006, Guidetti et al. 2008). A key question is whether the 
relation between magnetic field strength and density in galaxy 
groups is a continuation of this trend. 

This paper presents a detailed analysis of Faraday rotation 
in 3C 449, a bright, extended radio source hosted by the central 
galaxy of a nearby group. With the aim of shedding new light 
on the environment around this source, we derive the statistical 
properties of the magnetic field from observations of Faraday ro- 
tation, following the method developed by Murgia et al. (2004). 
We use numerical simulations to predict the Faraday rotation for 
dififerent strengths and power spectra of the magnetic field. 

The paper is organized as follows. In Sect.|2] the general 
properties of the radio source under investigation are presented. 
Sect.Opresents the radio images on which our analysis is based. 
In Sect.m we discuss the observed Faraday rotation distribu- 
tion of 3C 449 and assess the contribution from our Galaxy. The 
observed depolarization and its relation to the RM properties 
are investigated in Sect.|5] Our two-dimensional analysis of the 
structure of the RM fluctuations and a three-dimensional model 
of the magnetic field consistent with these results are presented 
in Sect. |6] and Sect.|7] respectively. Sect. |8] summarizes our con- 
clusions and briefly compares the Faraday -rotation properties of 
3C449 with those of other sources. 

Throughout this paper we assume a cosmology with //q = 
71 km s 'Mpc ', Q„, - 0.3, and Oa - 0.7, which implies that 
1 arcsec corresponds to 0.342kpc at the distance of 3C449. 



2. The radio source 3C449: general properties 

In this paper, we image and model the Faraday rotation distribu- 
tion across the giant Fanaroff'-Riley Class I (FRI; Fanaroff' & 
Riley 1974) radio source 3C449, whose environment is very 
similar to that of 3C31. The optical counterpart of 3C449, 



' ^ = 11/3 is the slope of the three-dimensional power spectrum for 
Kolmogorov turbulence. 



UGC 12064, is a dumb-bell galaxy and is the most promi- 
nent member of the group of galaxies 2231.2H-3732 (Zwicky 
& Kowal 1968). The source is relatively nearby (z=0. 017085, 
RC3.9, De Vaucouleurs et al. 1991) and quite extended, both in 
angular (30') and linear size, so it is an ideal target for an anal- 
ysis of the Faraday rotation distribution: detailed images can be 
constructed that can serve as the basis of an accurate study of 
magnetic field power spectra. 

3C 449 was one of the first radio galaxies studied in detail 
with the VLA (Perley et. al 1979). High- and low-resolution ra- 
dio data already exist and the source has been mapped at many 
frequencies. The radio emission of 3C449 (Fig.[T]) is elongated 
in the N-S direction and is characterized by long, two-sided jets 
with striking mirror symmetry close to the nucleus. The jets ter- 
minate in well-defined inner lobes, which fade into well polar- 
ized spurs, of which the southern one is more collimated. The 
spurs in turn expand to form diffuse outer lobes. 

The brightness ratio of the radio jets is very nearly 1, imply- 
ing that they are close to the plane of the sky if they are intrinsi- 
cally symmetrical and have relativistic flow velocities similar to 
those derived for other FRI jets (Perley et al. 1979, Feretti et al. 
1999a, Laing & Bridle 2002). In this paper, we therefore assume 
that the jets lie exactly in the plane of the sky, simplifying the 
geometry of the Faraday-rotating medium. 

Hot gas associated with the galaxy has been detected on both 
the group and galactic scales by X-ray imaging (Hardcastle et al. 
1998; Croston et al. 2003). These observations revealed deficits 
in the X-ray surface brightness at the positions of the outer ra- 
dio lobes, suggesting interactions with the surrounding material. 
Fig-fflshows radio contours at 1.365 GHz overlaid on the X-ray 
emission as observed by the XMM-Newton satellite (Croston et 
al. 2003). The X-ray radial surface brightness profile of 3C449 
derived from these data can be fitted with the sum of a point- 
source convolved with the instrumental response and a /3 model 
(Cavaliere & Fusco-Femiano, 1976), 

n,{r)^no{l+r^/r^y^'', (1) 

where r, and no are the distance from the group X-ray cen- 
tre, the group core radius and the central electron density, re- 
spectively. Croston et al. (2008) found a best fitting model with 
/5 = 0.42 + 0.05, r, = 57.1 arcsec and «o = 3.7 x lO^^cm^l In 
what follows, we assume that the group gas density is described 
by the model of Croston et al. (2008). The X-ray depressions 
noted by Croston et al. (2003) are at distances larger than those at 
which we can measure linear polarization, and there is no direct 
evidence of smaller cavities close to the nucleus. We therefore 
neglect any departures from the spherically-symmetrical density 
model, noting that this approximation may become increasingly 
inaccurate where the source widens (i.e. in the inner lobes and 
spurs). 

3C 449 resembles 3C 3 1 in environment and in radio mor- 
phology: both sources are associated with the central members 
of groups of galaxies and their redshifts are very similar. The 
nearest neighbours are at a projected distances of about 30 kpc in 
both cases. Both radio sources have large angular extents, bend- 
ing jets and long, narrow tails with low surface brightnesses and 
steep spectra, although 3C 3 1 appears much more distorted on 
large scales. There is one significant difference: the inner jets 
of 3C 3 1 are thought to be inclined by a;50° to the line of sight 
(Laing & Bridle 2002), whereas those in 3C449 are likely to 
be close to the plane of the sky (Feretti et al. 1999a). We there- 
fore expect that the magnetized foreground medium will be very 
similar in the two sources, but that the geometry will be signif- 
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Fig. 1. Radio contours of 3C449 at 1.365 GHz superposed on the XMM-Newton X-ray image (courtesy of J. Croston, Croston et 
al. 2003). The radio contours start at 3cr/ and increase by factors of 2. The restoring beam is 5.5 arcsec FWHM. The main regions 
of 3C 449 discussed in the text are labelled. 



icantly different, leading to a much more symmetrical distribu- 
tion of Faraday rotation in 3C 449 compared with that observed 
in 3C 31 by Laing et al. (2008). 

3. Total intensity and polarization properties 

The Very Large Array (VLA0) observations and their reduction 
have been presented by Feretti et al. (1999a). The high quality 
of these data make this source suited for a very detailed analysis 
of the statistics of the Faraday rotation. 

We produced total intensity (/) and polarization (Q and U) 
images at frequencies in the range 1.365 - 8.385 GHz from the 
combined, self-calibrated u-v datasets described by Feretti et al. 
(1999a). The centre frequencies and bandwidths are listed in 
Table [1] Each frequency channel was imaged separately, except 
for those at 8.245 and 8.445 GHz, which were averaged. The 
analysis below confirms that these frequency-bandwidth combi- 
nations lead to negligible Faraday rotation across the channels, 
as already noted by Feretti et al. (1999a). All of the datasets 
were imaged with Gaussian tapering in the u-v plane to give res- 
olutions of 1.25 arcsec and 5.5 arcsec FWHM, CLEANed and 
restored with circular Gaussian beams. The first angular res- 
olution is the highest possible at all frequencies and provides 
good signal-to-noise for the radio emission within 150 arcsec 
(^50kpc) of the radio core (the well defined radio jets and the 

^ The Very Large Array is a facility of the National Science 
Foundation, operated under cooperative agreement by Associated 
Universities, Inc. 



inner lobes), while minimizing beam depolarization. The lower 
resolution of 5.5 arcsec allows imaging of the extended emission 
as far as 300 arcsec (^100 kpc) from the core at frequencies from 
1.365 - 4.985 GHz (the 8.385-GHz dataset does not have ade- 
quate sensitivity to image the outer parts of the source). We can 
therefore study the structure of the magnetic field in the spur re- 
gions, which lie well outside the bulk of the X-ray emitting gas. 
Noise levels for both sets of images are given in Table [1] Note 
that the maximum scales of structure which can be imaged reli- 
ably with the VLA at 8.4 and 5 GHz are «180 and a!300 arcsec, 
respectively (Ulvestad, Perley & Chandler 2009). For this rea- 
son, we only use the Stokes / images for quantitative analysis 
within half these distances of the core. The Q and U images 
have much less structure on such large scales and are reliable to 
distances of +150 arcsec at 8.4 GHz and +300 arcsec at 5 GHz, 
limited by sensitivity rather than systematic errors due to miss- 
ing flux as in the case of / image. 

Images of polarized intensity P = {Q^ + U^y^^ (coiTected for 
Ricean bias, following Wardle & Kronberg 1974), fractional po- 
larization /:> = P// and polarization angle *F = (l/2)arctan(f//2) 
were derived from the /, Q, and U images. All of the polariza- 
tion images (P, p, ^') at a given frequency were blanked where 
the rms error in *F > 10° at any frequency. We then calculated 
the scalar mean degree of polarization (p) for each frequency 
and resolution; the results are listed in Table[T] The values of (p) 
are higher at 5.5 arcsec resolution than at 1.25 arcsec because of 
the contribution of the extended and highly polarized emission 
which is not seen at the higher resolution. At 1.25 arcsec, where 



4 



D. Guidetti et al.: Structure of the magnetoionic medium around the FR Class I radio galaxy 3C449 



Table 1. Parameters of the total intensity and polarization images. Col. 1: Observation frequency. Col. 2: Bandwidth (note that 
the images at 8.385 GHz are derived from the average of two frequency channels, both with bandwidths of 50 MHz, centred on 
8.285 and 8.485 GHz). Cols. 3, 4 : rms noise levels in total intensity (cr/) and linear polarization {ctqu, the average of ctq and ctu) 
at 1.25 arcsec FWHM resolution; Col. 5: mean degree of polarization at 1.25 arcsec; Col. 6, 7: rms noise levels for the 5.5 arcsec 
images; Col. 8: mean degree of polarization at 5.5 arcsec. We estimate that the uncertainty in the degree of polarization, which is 
dominated by systematic deconvolution errors on the / images, is a;0.02 at each frequency. 
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the beam depolarization is minimized, the mean fractional polar- 
ization shows a steady increase from 1.365 to 4.685 GHz, where 
it reaches an average value of 0.32 and then remains roughly 
constant at higher frequencies, suggesting that the depolariza- 
tion between 4.685 and 8.385 GHz is insignificant. 

4. The Faraday rotation in 3C449 

4.1. Rotation measure images 

A magnetized, ionized medium rotates the plane of polarization 
of linearly polarized radiation passing through it as follows: 

A>P = »P(/l) - ^Fo = RM ^^ (2) 

where *P(/l) is the position angle observed at a wavelength A and 
*Po is the intrinsic position angle. The rotation measure (RM) is 
related to the electron density (n^,), the magnetic field along the 
line-of-sight (By), and the path-length (L) through the Faraday- 
rotating medium according to: 

RM [ja(j/n,2] = 812 I rig [cm-^^Bi^ ij,c]dl . (3) 
Jo 

Images of rotation measure can be obtained for radio sources 
by fitting to the polarization angle as a function of A^, taking into 
account the well-known problem of n;r ambiguities in the ob- 
served *F, as is done for example by the AIPS task RM. Removal 
of these ambiguities requires observations at many wavelengths 
well-spaced in A^. 

We produced images of RM and its associated rms error with 
resolutions of 1.25 arcsec and 5.5 arcsec (Fig.|2^ and b) using 
a version of the AIPS task RM modified by G. B. Taylor The 
1.25 arcsec-RM map was made by combining the maps of polar- 
ization E-vector (^F) at all seven frequencies available to us, so 
our sampling of A^ is very good. The RM map was calculated us- 
ing a weighted least-squares fit at pixels with polarization angle 
uncertainties <10°at all frequencies. It is essentially the same 
as the RM image of Feretti et al. (1999a), but with more strin- 
gent blanking. The average fitting error is ^1.4radm"^ and is 
almost constant over the whole RM image. The image of RM at 
5.5 arcsec resolution was produced using the polarization posi- 
tion angles at the 6 frequencies between 1.365 and 4.985 GHz 
(see Table[T]), using the same blanking criterion as at higher res- 
olution. 



Patches with different size are apparent in the 1.25 arcsec 
resolution map, with fluctuations down to scales of a few kpc. 
The bulk of the RM values range from about -220radm"^ 
up to -90radm"^, dominated by the Galactic contribution (see 
Sect. I4!2b. The RM distribution peaks at -161.7 radm"^, with 
rms dispersion ctj^m =19.7 radm"^. Note that we have not cor- 
rected the values of o-rm for the fitting error crjiMa- A first order 
correction would be crRM„„,. - (o"rm ~ "^rMi Given the low 
value for ctrmh,, the effect of this correction would be very small. 

As was noted by Feretti et al. (1999a), the RM distribution 
in the inner jets is highly symmetric about the core with RM 
^ -197 radm"^ at distances $15 arcsec. The symmetry of the 
RM distribution in the jets is broken at larger distances from 
the core: while the RM structure in the southern jet is homoge- 
neous, with values around -130radm"^, fluctuations on scales 
of ^lOarcsec (^ 3 kpc) around a (RM) of -160radm"^ are 
present in the northern jet. In both lobes, we observe similar 
patchy RM structures with mean values (RM) ^ -164radm"^ 
and (Trm - 1 6 rad nT^ . 

At 5.5 arcsec resolution, more extended polarized regions of 
3C449 can be mapped with good sampling in A^. The aver- 
age fitting error is ^l.Oradm"^. Both the spurs are character- 
ized by (RM) ^ -160radm"^, with ctrm =15 and lOradm"^ 
in the north and south, respectively. The overall mean and 
rms for the 5.5 arcsec image, (RM) = -160.7 rad m"^ and 
CRM = 18.9 rad m"^, are very close to those determined at higher 
resolution and consistent with the integrated value of -162 + 
1 radm"^ derived by Simard-Normandin et al. (1981). 

It was demonstrated by Feretti et al. (1999a) that the polar- 
ization position angles at 1 .25 arcsec resolution accurately fol- 
low the relation A^F oc A^ over a large range of rotation. We find 
the same effect at lower resolution: plots of E-vector position 
angle 'I' against A^ at representative points of the 5.5 arcsec-RM 
image are shown in Fig[3] As at the higher resolution, there are 
no significant deviations from the relation A^P oc A^ over a range 
of rotation A*!* of 600 °, confirming that a foreground magnetized 
medium is responsible for the majority of the Faraday rotation 
and extending this result to regions of lower surface brightness. 

In Fig. [2c) and (d), we show profiles of ctrm for both low 
and high resolution RM images. The 1 .25 arcsec profile was ob- 
tained by averaging over boxes with lengths ranging from 9 to 
13 kpc along the radio axis; for the 5.5 arcsec profile we used 
boxes with a fixed length of 9 kpc (these sizes were chosen to 
give an adequate number of independent points per box). The 
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Fig. 2. (a): Image of the rotation measure of 3C449 at a resolution of 1.25arcsec FWHM, computed at the seven frequencies 
between 1.365 and 8.385 GHz. (b): Image of the rotation measure of 3C449 at a resolution of 5.5 arcsec FWHM, computed at the 
six frequencies between 1.365 and 4.985 GHz. In both of the RM images, the sub-regions used for the two-dimensional analysis 
of Sect.|6]are labelled, (c) and (d): profiles of ctrm as a function of the projected distance from the radio source centre. The points 
represent the values of ctrm evaluated in boxes as described in the text. The horizontal and vertical bars represent the bin widths and 
the rms on the mean expected from fitting errors, respectively. Positive distances are in the direction of the north jet and the vertical 
dashed lines show the position of the nucleus. 



boxes extend far enough perpendicular to the source axis to in- 
clude all unblanked pixels. In both plots, there is clear evidence 
for a decrease in the observed ctrm towards the periphery of the 
source, the value dropping from ^30rad nT^ close to the nucleus 
to ^lOradm"^ at 50kpc. This is qualitatively as expected for 
foreground Faraday rotation by a medium whose density (and 
presumably also magnetic field strength) decreases with radius. 
The symmetry of the o-f^M profiles is consistent with our assump- 
tion that the radio source lies in the plane of the sky. 



4.2. The Galactic Faraday rotation 

For the purpose of this work, 3C449 has an unfortunate line 
of sight within our Galaxy. Firstly, the source is located at I - 
95.4° = -15.9° in Galactic coordinates, where the Galactic 
magnetic field is known to be aligned almost along the line of 
sight. Secondly, there is evidence from radio and optical imaging 
for a diffuse, ionized Galactic feature in front of 3C 449, perhaps 
associated with the nearby HII region S126 (Andernach et al. 
1992). Estimates of the Galactic foreground RM at the position 
of 3C 449 from observations of other radio sources are uncertain: 
Andernach et al. (1992) found a mean value of -212radm"^ for 
six nearby sources, but the spherical harmonic models of Dineen 



& Coles (2005), which are derived by fitting to the RM values 
of large numbers of extragalactic sources, predict -135radm"^. 
Nevertheless, it is clear that the bulk of the mean RM of 3C 449 
must be Galactic. 

In order to investigate the magnetized plasma local to 
3C 449, we need to constrain the value and possible spatial vari- 
ation of this Galactic contribution. The profiles of ctrm (Fig. |2]i 
show that the small-scale fluctuations of RM drop rapidly 
with distance from the nucleus. We might therefore expect the 
Galactic contribution to dominate on the largest scales. At low 
resolution, we can determine the RM accurately out to wlOOkpc 
from the core. This is roughly 5 core radii for the X-ray emission 
and therefore well outside the bulk of the intra-group gas. 

In order to estimate the Galactic RM contribution, we aver- 
aged the 5.5-arcsec RM image in boxes of length 20kpc along 
the radio axis (the box size has been increased from that of Fig.|2] 
to improve the display of large-scale variations). The profile of 
(RM) against the distance from the radio core is shown in Fig.|4l 
The large deviations from the mean in the innermost two bins 
are associated with the maximum in ctrm and are almost cer- 
tainly due to the intra-group medium. The dispersion in (RM) is 
quite small in the south and the value of (RM) = -160.7radm"^ 
for the whole source is very close to that of the outer south jet. 
There are significant fluctuations in the north, however. Given 
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Fig. 3. Plots of E-vector position angle ^' against at repre- 
sentative points of the 5.5-arcsec RM map. Fits to the relation 
= *Po + RM/l^ are shown. The values of RM are given in 
the individual panels. 



their rather small scale (~300 arcsec), it is most likely that these 
arise in the local environment of 3C 449 and we include them in 
the statistical analysis given below. 

There is some evidence for linear gradients in Galactic RM 
on arcminute scales: Laing et al. (2006) found a gradient of mag- 
nitude 0.025 radm"^ arcsec"' along the jets of the radio galaxy 
NGC 315 (/ = 124.6 ° -32.5 °). They argued that this gra- 
dient is almost certainly Galactic in origin, since the amplitude 
of the linear variation exceeds that of the small-scale fluctua- 
tions associated with NGC 315. In order to check the eff'ect of 
a large-scale Galactic RM gradient on our results, we computed 
an unweighted least-squares fit of a function (RM) - RMq + ax, 
where a and RMq are constant and x is measured along the 
radio axis. The two innermost bins in Fig.|4] were excluded 
from the fit. Our best estimate for the gradient is very small: 
a=0.0054radm"^ arcsec"'. We have verified that subtraction of 
this gradient has a negligible effect on the structure-function 
analysis given in Sect. 17.31 

We therefore adopt a constant value of - 160.7 rad nT^ as the 
Galactic contribution. 



5. Depolarization 

Faraday rotation generally leads to a decrease of the degree of 
polarization with increasing wavelength, or depolarization. We 
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The black vertical dashed line indicates the position of the nu- 
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define DP^^ = p(Ai)/ p(A2), where p(A) is the degree of polariza- 
tion at a given wavelength A. We adopt the conventional usage 
in which higher depolarization corresponds to a lower value of 
DP 

Laing (1984) has summarized the interpretation of polariza- 
tion data. Faraday depolarization of radio emission from radio 
sources can occur in three principal ways: 

1 . thermal plasma is mixed with the synchrotron emitting ma- 
terial (internal depolarization); 

2. there are fluctuations of the foreground Faraday rotation 
across the beam (beam depolarization) and 

3. the polarization angle varies across the finite band of the re- 
ceiving system (bandwidth depolarization). 

We first estimated the bandwidth effects on the polar- 
ized emission of 3C449 using the RM measurements from 
Sect. 14.11 In the worst case (the highest absolute RM value of 
-240 rad m"^) at the lowest frequency of 1.365 GHz) the rota- 
tion across the band is ^10°. This results in a depolarization of 
1.7%, negligible compared with errors due to noise. 

If A^ rotation is observed over a position-angle range »90°, 
then a foreground screen must be responsible for the bulk of the 
observed RM. In that case, depolarization can still result from 
unresolved inhomogeneities of thermal density or magnetic field 
in the surrounding medium. Our analysis of the depolarization 
of 3C449 is based on the approach of Laing et al. (2008). In the 
presence of a foreground Faraday screen with a small gradient 
of RM across the beam, it is still possible to observe A^ rotation 
over a wide range of polarization angle and the wavelength de- 
pendence of the depolarization is expected to follow the Burn 
law (Burn 1966): 



p(A) = p(Q)cxp(-kA\ 



(4) 
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where /:i(0) is the intrinsic value of the degree of polarization and 
k=2\VRM\W, with FWHM = 2cr(2 ln2)i/2. Since k oc \VRM\^, 
Eq.|4]clearly illustrates that higher RM gradients across the beam 
generate higher k values and hence higher depolarization. The 
variation of p with wavelength can potentially be used to es- 
timate fluctuations of RM across the beam which are below the 
resolution limit. We can determine the intrinsic polarization p{0) 
and the proportionality constant A; by a linear fit to the logarithm 
of the observed fractional polarization as a function of A'^. 

We made images of k at both standard resolutions by 
weighted least-squares fitting to the fractional polarization maps, 
using the FARADAY code by M. Murgia. The same frequencies 
were used as for the RM images: 8.385 - 1.365 GHz and 4.985 
- 1.365 GHz at 1.25 and 5.5 arcsec resolution, respectively. By 
simulating the error distributions for p, we established that the 
mean values of k were biased significantly at low signal-to-noise 
(cf. Laing et al. 2008), so only data with p > AcTp at each fre- 
quency are included in the fits. We estimate that any bias is neg- 
ligible compared with the fitting error. We also derived profiles 
of k using the same sets of boxes as for the ctrm profiles in Fig.|2] 

The 1.25 arcsec resolution k map is shown in Fig.|5la), to- 
gether with the profile of the k values (Fig.|5};). The fit to a law 
is very good everywhere: examples of fits at selected pixels in the 
jets and lobes are shown in in Fig.|6] The symmetry observed in 



the ctrm profiles is also seen in the 1.25 arcsec k image (Fig.|5]l: 
the mean values of k are ^50rad^m"'* for both lobes, 107 and 
82 rad^ m for the northern and southern jet, respectively. The 
region with the highest depolarization is in the northern jet, very 
close to the core and along the west side. The integrated value 
of k at this resolution is ^56 rad^ m""*, corresponding to a mean 
depolarization DP^l^ ^ 0.87. 

The image and profile of k at 5.5-arcsec resolution are shown 
in Fig.|5lb) and (d). The fit to a A^ law is in general good and ex- 
amples are shown in in Figs.|6] As mentioned earlier, the max- 
imum scale of structure imaged accurately in total intensity at 
5 GHz is ~300 arcsec (100 kpc) and there are likely to be sig- 
nificant systematic errors in the degree of polarization on larger 
scales. We therefore show the profile only for the inner +50 kpc. 
Over this range, the k profiles are quite symmetrical, as at higher 
resolution. Note also that the small regions of very high k at the 
edge of the northern and southern spurs in the map shown in 
Fig.|5]are likely to be spurious. 

The mean values of k are ^184rad^m"^ and 178rad^m"^ 
for the northern and southern lobes, respectively; and ^238 and 
174 rad^ m "* in the northern and in the southern spurs. The inte- 
grated value of k is ^194rad^m corresponding to a depolar- 
ization DP20™ ~ 0.64. 
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Fig. 7. Plots of degree of polarization, p (log scale) against /l"* 
for representative points at 5.5-arcsec resolution. Burn law fits 
(equation© are also plotted. The values of k are quoted in the 
individual panels. 



To summarize, we observe depolarization between 20 cm 
and 3 cm. Since we measure lower values of ^ at 1 .25 arcsec than 
5.5 arcsec, there is less depolarization at high resolution, as ex- 
pected in the case of beam depolarization. The highest depolar- 
ization is observed in a region of the northern jet, close to the ra- 
dio core and associated with a large RM gradient. Depolarization 
is significantly higher close to the nucleus, consistent with the 
higher path length through the group gas observed in X-rays. 



Aside from this global variation, we have found no evidence 
for detailed correlation of depolarization with source structure. 
Depolarization and RM data are therefore both consistent with a 
foreground Faraday screen. We show in Sect. l6.2l that the resid- 
ual depolarization at 1.25-arcsec resolution can be produced 
by RM fluctuations on scales smaller than the beamwidth, but 
higher-resolution observations are needed to establish this con- 
clusively. 

6. Two Dimensional Analysis 

6.1. General considerations 

In order to interpret the fluctuations of the magnetic field re- 
sponsible for the observed RM and depolarization of 3C449, 
we first discuss the statistics of the RM fluctuations in two di- 
mensions. We use the notation of Laing et al. (2008), in which 
f = (/i,/v,/z) is a vector in the spatial frequency domain, corre- 
sponding to the position vector r - (x,y,z). We take the z-axis 
to be along the line of sight, so that the vector = (x, y) is in 
the plane of the sky and f - {fx, fy) is the corresponding spatial 
frequency vector. Our goal is to estimate the RM power spec- 
trum C(fj.), where C(f±)dfxdfy is the power in the area df-dfy 
and in turn to derive the three dimensional magnetic field power 
spectrum w(f), defined so that w(f)dfxdfydf, is the power in a 
volume dfxdfydf; of frequency space. 

The relation between the magnetic field statistics and the ob- 
served RM distribution is in general quite complicated, depend- 
ing on fluctuations in the thermal gas density, the geometry of 
the source and the surrounding medium and the effects of in- 
complete sampling. In order to derive the magnetic-field power 
spectrum, we make the following simplifying assumptions, as in 
Guidetti et al. (2008) and Laing et al. (2008). 

1 . The observed Faraday rotation is due entirely to a foreground 
ionized medium (in agreement with our results in Sects|4]and 

E). 

2. The magnetic field is an isotropic, Gaussian random variable, 
and can therefore be characterized by a power spectrum w{f) 
which is a function only of scalar frequency /. 

3. The form of the magnetic field power spectrum is indepen- 
dent of position. 

4. The magnetic field is distributed throughout the Faraday- 
rotating medium, whose density is a smooth, spherically 
symmetric function. 

5. The amplitude of w{ f) is spatially variable, but is a function 
only of the thermal electron density. 

These assumptions guarantee that the spatial distribution of the 
magnetic field can be described entirely by its power spectrum 
w{ f), and that for a medium of constant depth and density, the 
power spectra of magnetic field and RM are proportional (EnBlin 
and Vogt 2003). 

If the fluctuations are isotropic, the RM power spectrum 
C(/j.) is the Hankel transform of the autocorrelation function 
C{r^^, defined as 



C(r^) = (RM(r^ + rl)RM(rl)), 



(5) 



where and r'j^ are vectors in the plane of the sky and () is an 
average over r^. 

In an ideal case, it would be possible to derive the RM power 
spectrum and, consequently, that of the magnetic field, directly 
from C(rx). In reality, the observations are affected first by the 
effects of convolution with the beam, which modify the spatial 
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Fig. 8. (a)-(f): Plots of the RM structure functions for the sub-regions showed in Fig|2] The horizontal bars represent the bin widths 
and the crosses the centroids for data included in the bins. The red lines are the predictions for the CPL power spectra described 
in the text, including the effects of the convolving beam. The vertical eiTor bars are the rms variations for the structure functions 
derived using a CPL power spectrum with the quoted value of q on the observed grid of points for each sub-region, (g)-(l): as (a)-(f) 
but using a BPL power spectra with fixed slopes and break frequency, but variable normalization. 



statistics of RM, and secondly, by the limited size and irregu- 
lar shape of the sampling region for 3C 449, which results in a 
complicated window function (EnBlin & Vogt, 2003) and lim- 
its the accuracy with which the zero-level can be determined. 
In Sect. 14. 21 we showed that the Galactic contribution to the 
3C449 RM is large, and we argued that a constant value of 
-160.7radm"^ is the best estimate for its value. Fluctuations in 
the Galactic magnetic field on scales comparable with the size 
of the radio sources could be present; conversely, the local envi- 
ronment of the source might make a significant contribution to 
the mean RM. Both of these possibilities lead to difficulties in 
the use of the autocoiTelation function. 

Laing et al. (2008) demonstrated a procedure which takes 
into account the convolution effects and minimises the effects 
of uncertainties in the zero-level. In particular, they showed that: 

1. In the short-wavelength limit (meaning that changes in 
Faraday rotation across the beam are adequately represented 
as a linear gradient), the measured RM distribution is closely 
approximated by the convolution of the true RM distribution 
with the observing beam. 



2. The structure function is a powerful and reliable statistical 
tool to quantify the two dimensional fluctuations of RM, 
given that it is independent of the zero level and structure 
on scales larger than the area under investigation. 

The structure function is defined by 

S{r^) = < [RM(r^ + r^) - RM(r'J]^ > (6) 

(Simonetti, Cordes & Spangler, 1984; Minter and Spangler 
1996). It is related to the autocoiTelation function C{rj_) for a 
sufficiently large averaging region by S (r^) - 2[C(rx) - C(0)]. 

Laing et al. (2008) also derived the effects of convolution 
with the observing beam on the observed structure function. 
For the special case of a power-law power spectrum (their Eq. 
B2), they showed that the observed structure function after con- 
volution can be heavily modified even at separations up to a 
few times the FWHM of the observing beam. This effect must 
be taken into account when comparing observed and predicted 
structure functions. Laing et al. (2008) and Guidetti et al. (2008) 
also showed that numerical simulation of depolarization pro- 
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vides complementary information on RM fluctuations on scales 
smaller than the beam. 

Following the approach of Laing et al. 2008, we initially 
used the RM structure function to determine the form of w(f) 
(Sect. I6.2I 1. while for its normalization (determined by global 
variations of density and magnetic field strength), we made use 
of three-dimensional simulations (Sect.|7]i. 



6.2. Structure functions 

We calculated the structure function for discrete regions of 
3C449 over which we expect the spatial variations of thermal 
gas density, rms magnetic field strength and path length to be 
reasonably small. For each of these regions, we first made un- 
weighted fits of model structure functions derived from power 
spectra with simple, parameterized functional forms, account- 
ing for convolution with the observing beam. We then generated 
multiple realizations of a Gaussian, isotropic, random RM field, 
with the best-fitting power spectrum on the observed grids, again 
taking into account the effects of the convolving beam. Finally, 
we made a weighted fit using the dispersion of the synthetic 
structure functions as estimates of the statistical errors for the 
observed structure functions, which are impossible to quantify 
analytically (Laing et al. 2008). These errors, which result from 
incomplete sampling, are much larger than those due to noise, 
but depend only weakly on the precise form of the underlying 
power spectrum. Our measure of goodness of fit is x^, summed 
over a range of separations from rj_ = FWHM to roughly half 
of the size of the region: there is no information in the structure 
function for scales smaller than the beam, and the upper limit is 
set by sampling. The errors are, of course, much higher at the 
large spatial scales, which are less well sampled. Note, however, 
that estimates of the structure function from neighbouring bins 
are not statistically independent, so it is not straightforward to 
define the effective number of degrees of freedom. 

We selected six regions for the structure-function analysis, 
as shown in Fig. [2] These are symmetrically placed about the 
nucleus, consistent with the orientation of the radio jets close 
to the plane of the sky. For the north and south jets, we de- 
rived the structure functions only at 1 .25-arcsec resolution, as the 
low-resolution RM image shows no additional structure and has 
poorer sampling. For the north and south lobes, we computed 
the structure functions at both resolutions over identical areas 
and compared them. The agreement is very good, and the low- 
resolution RM images do not sample significantly larger spatial 
scales, so we show only the 1. 25-arcsec results. Finally, we used 
the 5.5-arcsec RM images to compute the structure functions for 
the north and south spurs, which are not detected at the higher 
resolution. 

The structure function has a positive bias given by 2cr^^.^^ 
where <T„oise is the uncorrected random noise in the RM im- 
age (Simonetti, Cordes& Spangler, 1984). The mean noise of the 
1.25 and 5.5-arcsec RM maps is <1 radm^ and essentially uncor- 
rected on scales larger than the beam. For each region we there- 
fore subtracted 2cr^ . from the structure functions, although this 

noise ' ^ 

correction is always small. The noise-corrected structure func- 
tions are shown in Fig. [8] 

The individual observed structure functions have approxi- 
mately power-law forms. Given that the structure function for 
a power-law power spectrum with no frequency limits is itself a 
power law (Minter & Spangler 1996; Laing et al. 2008), we first 



tried to fit the observed data with a RM power spectrum of the 
form 

C(A) <^ A'' (7) 

over an infinite frequency range. This last assumption allows us 
to use the analytical solution of the structure function, including 
convolution (Laing et al. 2008) and therefore to avoid numerical 
integration. 

The fits were quite good, but systematically gave slightly too 
much power on small spatial scales and over-predicted the de- 
polarization. We therefore fit a cut-off power law (CPL) power 
spectrum 

C(A) = A < 

- Cof^ f± < /max 

= A > /max . (8) 

Initially, we consider values of /min sufficiently small that their 
effects on the structure functions over the observed range of 
separations are negligible. The free parameters of the fit in this 
case are the slope, q, the cut-off spatial frequency /^ax and the 
normalization of the power spectrum, Cq. In Table|2] we give 
the best-fitting parameters for CPL fits to all of the individ- 
ual regions. The fitted model structure functions are plotted in 
Fig.[8ja)-(f), together with error bars derived from multiple re- 
alizations of the power spectrum as in Laing et al. (2008). 

In order to constrain RM structure on spatial scales below 
the beamwidth, we estimated the depolarization expected from 
the best power spectrum for each of the regions with 1 .25-arcsec 
RM images, following the approach of Laing et al. (2008). To do 
this, we made multiple realizations of RM images on an 8192^ 
grid with fine spatial sampling. We then derived the Q and U 
images at our observing frequencies, convolved to the appropri- 
ate resolution and compared the predicted and observed mean 
degrees of polarization. These values are given in Table |4] The 
uncertainties in the expected < k > in Table |4]represent statisti- 
cal errors determined from multiple realizations of RM images 
with the same set of power spectrum parameters. The predicted 
and observed values are in excellent agreement. A constant value 
of /max - 1.67arcsec"' predicts very similar values, also listed 
in Table |4] We have not compared the depolarization data at 
5.5-arcsec resolution in the spurs because of limited coverage 
of large spatial scales in the / images (Sect.|3]l, which is likely to 
introduce systematic errors at 4.6 and 5.0GHz. 

We performed a joint fit of the CPL power spectra, minimiz- 
ing the summed over all six sub-regions, giving equal weight 
to each and allowing the normalizations to vary independently. 
In this case the free parameter of the fit are: the six normal- 
izations (one for each sub-region) the slope and the maximum 
spatial frequency. The joint best-fitting single power-law power 
spectrum has q = 2.68. 

A single power law slope does not give a good fit to all of 
the regions simultaneously, however. It is clear from Fig. [8] and 
Table|2]that there is a flattening in the slope of the observed struc- 
ture functions on the largest scales (which are sampled primarily 
by the spurs). In order to fit all of the data accurately with a sin- 
gle functional form for the power spectrum, we adopt a broken 
power law form (BPL) for the RM power spectrum: 

C(A) = A < /m,„ 

= Dof^"'-""'/-'' ft > A 

= Dofj_ /max > f± > fh 

= A > /max . (9) 
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Table 2. CPL power spectrum parameters for the six individual sub-regions of 3C 449 (lower and upper limits are quoted at -90% 
confidence). 
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We performed a BPL joint fit, in the same way as for the CPL 
power spectra. In this case the free parameters of the fit are: the 
six normalizations. Do, one for each sub-region, the high and 
low-frequency slopes, q\i and qi, the break and maximum spatial 
frequency fi, and /max- We found best fitting parameters of qi- 
2.07, ^h=2.98, /b=0.031 arcsec"'. As noted earlier, we also fixed 
/max - 1 -67 arcsec"' to ensure consistency with the observed de- 
polarizations at 1.25-arcsec resolution. The corresponding struc- 
ture functions are plotted in Fig.[8jg)-(1) and the normalizations 
for the individual regions are given in Table |4] As for the CPL 
fits, the errors bars are derived from the rms scatter of the struc- 
ture functions of multiple convolved RM realizations. 

It is evident from Fig.|8] that the structure functions corre- 
sponding to the BPL power spectrum, which gives less power 
on large spatial scales, are in much better agreement with the 
data. The joint BPL fit has a of 17.7, compared with 33.5 
for the joint CPL fit (the former has only two extra parameters) 
confirming this result. 

We have so far ignored the effects of any outer scale of 
magnetic -field fluctuations. This is justified by the fact that the 
structure functions for the spurs continue to rise at the largest 
observed separations, indicating that the outer scale must be 
^10 arcsec (^30kpc). The model structure functions fit to the 
observations assume that the outer scale is infinite and the real- 
izations are generated on sufficiently large grids in Fourier space 
that the effects of the implicit outer scale are negligible over the 
range of scales we sample. We use structure-function data for 
the entire source to determine an approximate value for the outer 
scale in Sect. 17.31 

We now adopt the BPL power spectrum with these parame- 
ters and investigate the spatial variations of the RM fluctuation 
amplitude using three-dimensional simulations. 

7. Three-dimensional analysis 

7.1. Models 

We used the software package FARADAY (Murgia et al. 
2004) to compare the observed RM with simulated images de- 
rived from three-dimensional multi-scale magnetic-field models. 
Given a field model and the density distribution of the thermal 
gas, FARADAY calculates an RM image by integrating Eq.[3] 
numerically. As in Sect. |6l we model the fluctuations of RM 
on the assumption that the magnetic field responsible for the 
foreground rotation is an isotropic, Gaussian random variable 
and therefore characterized entirely by its power spectrum. Each 
point in a cube in Fourier space is first assigned components of 
the magnetic vector potential. The amplitudes are selected from 
a Rayleigh distribution of unit variance and the phases are ran- 
dom in [0, 27r]. The amplitudes are then multiplied by the square 
root of the power spectrum of the vector potential, which is sim- 



ply related to that of the magnetic field. The coiTesponding com- 
ponents of the magnetic field along the line of sight are then cal- 
culated and transformed to real space. This procedure ensures 
that the magnetic field is divergence-free. The field components 
in real space are then multiplied by the model density distribu- 
tion and integrated along the line of sight to give a synthetic RM 
image at the full resolution of the simulation, which is then con- 
volved to the observing resolution. 

For 3C 449, we assumed that the source is in a plane perpen- 
dicular to the line of sight which passes through the group centre 
and simulated the field and density structure using a 2048^ cube 
with a real-space pixel size of 0.1 kpc. We used the best-fitting 
BPL power spectrum found in Sect. l6.2l but with a spatially- 
variable normalization, as described below (Sect. 17.2b . and a 
low-frequency cut-off /min, corresponding to a maximum scale 
of the magnetic field fluctuationsjj Amax (= /mil,)- The power 
spectrum of Eq. |9]is then set to for / < /min ■ We fixed the min- 
imum scale of the fluctuations Amin = 0.2 kpc. This is equivalent 
to the value /max = L67 arcsec"' found in Sect. l6.2l and also con- 
sistent with the requirement that the minimum scale can be no 
larger than twice the pixel size for adequate sampling. 

We made multiple synthetic RM images at resolutions of 
1.25 and 5.5 arcsec over the fields of view of the observations 
for each combination of parameters. In order to estimate the spa- 
tial variation of the magnetic-field strength, we first made a set 
of simulations with a large, fixed value of Amax and compared 
the predicted and observed profiles of ctrm (Sect. I7.2l i. We then 
fixed the radial variation of the field at its best-fitting form and 
estimated the value of Amax using a structure-function analysis 
for the whole source (Sect. l7.3l l. 

7.2. Magnetic field strength and radial profile 

In order to estimate the radial variation of field strength, we 
first fixed the value of the outer scale to be Amax - 205 kpc, 
the largest allowed by our simulation grid. Our approach was 
to make a large number of simulations for each combination of 
field strength and radial profile and to compare the predicted 
and observed values of ctrm evaluated over the boxes used in 
Sect. l4.1l (Fig.l2li. We used;^'^ summed over the boxes as a mea- 
sure of the goodness of fit. This procedure is independent of the 
precise value of the outer scale provided that it is much larger 
than the averaging boxes. We express our results in terms ofx^^^, 
which is the value of divided by the number of degrees of 
freedom. 



^ Here we refer to the scale length A as a complete wavelength, i.e. 
A = 1 //. This differs by a factor of 2 from the definition in Guidetti 
et al. 2008, where A is the reversal scale of the magnetic field, so A = 
1/2/. 



12 



D. Guidetti et al.: Structure of the magnetoionic medium around the FR Class I radio galaxy 3C449 



Table 3. Best-fitting parameters for the joint CPL and BPL fits to all six sub-regions of 3C449 (lower and upper limits are quoted 
at ~90% confidence). The values of q and /max for the joint CPL fit and qh, qi and /b for the joint BPL fit are the same for all 
sub-regions, while the normalizations are varied to minimize the overall In the joint BPL fit, the maximum frequency is fixed at 
/max - 1.67 arcsec"'. 







Best Fit 






Min Slope 




Max Slope 




q 


^max 




q 


fm<iX 


q" 


/nax 


joint CPL 


2.68 


1.67 


33.5 


2.55 


1.30 


2.81 


2.00 






Best Fit 






Min Slope 




Max Slope 




H 


/b qi, 




qj 


/b qJ, 


q^ 


/b q^ 


joint BPL 


2.07 


0.031 2.97 




1.99 


0.044 2.91 


2.17 


0.021 3.09 



Table 4. Normalizations Co and Do for the individual fit parameters corresponding to the CPL, joint CPL and joint BPL fits at a 
resolution of 1.25 arcsec. Observed and expected depolarizations are also given. Col.l: region; Col. 2: observed Burn law < k >. 
Col. 3, 4 and 5: normalization constant Co, fitted maximum spatial frequency /max for the best CPL power spectrum of each region 
and the predicted < k> for each power spectrum. Col. 6, 7 as Col. 3 and 5 but for the joint fit to each CPL power spectrum; Col. 8 
and 9 as Col. 3, and 5 but for the joint BPL power spectrum. For both the joint CPL and BPL fits, the maximum frequency is fixed 
at /max - L67 arcsec"'. In calculating each value of < A: > only data with p > AcTp are included. 



Region 


Observed <k> 




CPL 




JOINT CPL 


JOINT BPL 






Co 


fmdy. 


<k> 


Co 


<k> 


D(, 


<k> 




(rad^ m""*) 




(arcsec"') 


(rad^ m"'') 




(rad^ m-'') 




(rad^ m"*) 


NLOBE 


61±6 


0.96 


1.60 


63±3 


1.52 


66±3 


1.91 


52±4 


N JET 


106±12 


1.34 


1.21 


106±5 


4.76 


110+5 


0.5 


109±5 


S JET 


91±11 


1.50 


1.95 


70+5 


1.94 


73±5 


1.52 


65±4 


SLOBE 


50±5 


1.18 


1.68 


53+2 


1.28 


50±2 


2.20 


45±3 



Table 5. Summary of magnetic field power spectrum and density scaUng parameters. 



BPL power spectrum 


ft = 


qi =2.07 
^h=2.98 
0.031 arcsec"' 


low-frequency slope 
high-frequency slope 
break frequency (A/^ = l//i,=l 1 kpc) 




= 1.67 arcsec"' 
/,„„ fitted 


maximum frequency (An,in = l//n,ax=0.2kpc) 
minimum frequency (An,ax = l/fmm) 


Scaling of the magnetic field 




Bo fitted 


Average magnetic field at group centre 




77 fitted 


Magnetic field exponent of the radial profile: <S)(r) = S, 1860 ['^J 



We initially tried a radial field-strength variation of the form: 



1/2 



riQ 



(10) 



as used by Guidetti et al. (2008) and Laing et al. (2008). Here, 
Bq is the rms magnetic field strength at the group centre and 
ng(r) is the thermal electron gas density, assumed to follow the 
jS-model profile derived by Croston et al. (2008; see Sect.|2l). This 
functional form is consistent with other observations, analytical 
models and numerical simulations. In particular, rj = 2/3 corre- 
sponds to flux-freezing and r] - 1/2 to equipartition between 
thermal and magnetic energy. Dolag et al. (2001) and Dolag 
(2006) found rj ^ 1 from the correlation between the observed 
rms RM and X-ray surface brightness in galaxy groups and clus- 
ters and showed that this is consistent with the results of MHD 
simulations. 



We produced simulated RM images for each combination 
of Bo and 77 in the ranges 0.5 - 10//G in steps of 0.1 yuG and 
- 2 in steps of 0.01, respectively. We then derived the synthetic 
CRM profiles and, by comparing them with the observed one, cal- 
culated the unweighted We repeated this procedure 35 times 
at each angular resolution, noting the (Bq, rf) pair which gave the 
lowest in each case. These values are plotted in Fig.|9] As in 
earlier work (Murgia et al. 2004; Guidetti et al. 2008; Laing et 
al. 2008), we found a degeneracy between values of Bq and 77, in 
the sense that the fitted values are positively correlated, but there 
are clear minima in;^'^ at both resolutions. We therefore adopted 
the mean values of Bo and 77, weighted by 1 Ix^, as the best over- 
all estimates. These are also plotted in Fig.|9] as blue crosses. 
Although the central magnetic field strengths derived for the two 
RM images are consistent at the Icr level (Bo=2.8+0.5/vG and 
Bq- 4.1 +1.2yuG at 5.5 and 1.25-arcsec resolution, respectively). 
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the values of rj are not. The best-fitting values are 77=0. 0+0. 1 at 
5.5 arcsec FWHM and 77=0.8+0.4 at 1.25 arcsec FWHM. 

We next produced 35 RM simulations at each angular reso- 
lution by fixing Bq and rj at their best values for that resolution. 
This allowed us to calculate weighted ;^'^'s for the ctrm profiles, 
evaluating the errors for each box by summing in quadrature the 
rms due to sampling (determined from the dispersion in the re- 
alizations) and the fitting eiTor of the observations. These values 
are listed in Table|6] The observed and best-fitting model profiles 
at both the angular resolutions are shown in Fig.fTOl 

A model with 77 ^ at all radii is unlikely a priori: pre- 
vious work has found values of 0.5 ^ 77 $ 1 in other sources 
(Dolag 2006; Guidetti et al. 2008; Laing et al. 2008). The most 
likely explanation for the low value of 77 inferred from the low- 
resolution image is that the electron density distribution is not 
well represented by a y6-model (which describes a spherical and 
smooth distribution) at large radii. In support of this idea, Fig.[T] 
shows that the morphology of the X-ray emission is not spher- 
ical at large radii, but quite irregular. Croston et al. (2003) and 
Croston et al. (2008) pointed out that the quality of the fit of a 
single y6-model to the X-ray surface brightness profile was poor 
in the outer regions, suggesting small-scale deviations in the gas 
distribution. The single ^S-model gave a better fit to the inner re- 
gion of the X-ray surface brightness profile, where the polarized 
emission of 3C449 can be observed at 1.25-arcsec resolution. 
Our aim is to fit the cr-gM profiles at both resolutions with the 
same distribution of nf{r){B{r)^yi^ . In the rest of this subsection 
we assume that the density profile iieif) is still represented by 
the single y6-model, even though we have argued that it might 
not be appropriate in the outer regions of the hot gas distribu- 
tion. Although the resulting estimates of field strength at large 
radii may be unreliable, the fit is still necessary for the calcula- 
tion of outer scale described in Sec. 17. 31 which depends only on 
the combined spatial variation of density and field strength. 

The best-fitting model at 1.25-arcsec resolution, which is 
characterized by a more physically reliable 77, gives a very bad fit 
to the low resolution profile at almost all distances from the core 
(Fig.fTOh). Conversely, the model determined at 5.5-arcsec res- 
olution gives a very poor fit to the sharp peak in cr-gM observed 
within 20kpc of the nucleus at 1.25-arcsec resolution, where the 
radio and X-ray data give the strongest constraints (Fig.fTOb). We 
have also verified that no single intermediate value of 77 gives an 
adequate fit to the cr^u profiles at all distances from the nucleus. 

A better description of the observed cr-gM profile is provided 
by the empirical function: 



= Bo 



no 



no 



r> r^. 



(11) 



where 77int, 77out are the inner and outer scaling index of the mag- 
netic field and is the break radius. 

We fixed 77int=1.0 and ?7out=0.0, consistent with our initial 
results, in order to reproduce both the inner sharp peak and the 
outer flat decline of the ct-bm observed at the two resolutions, 
keeping as a free parameter. We made three sets of three- 
dimensional simulations for values of the outer scale Anjax= 205, 
65 and 20kpc. Anticipating the result of Sect. 17.31 we plot the 
results only for Amax= 65kpc, but the derived ctrm profiles are 
in any case almost independent of the value of the outer scale in 
this range. The new simulations were made only at a resolution 



of 5.5 arcsec, since the larger field of view at this resolution is 
essential to define the change in slope of the profile. 

In order to determine the best-fitting break radius, in 
Eq.[TT] we produced 35 sets of synthetic RM images for a grid of 
values of Bo and for each outer scale, noting the pair of val- 
ues which gave the minimum unweighted for the o-rm profile 
for each set of simulations. These values are plotted in Fig. [TT] 
which shows that there is a degeneracy between the break radius 
Tn, and Bo- As with the similar degeneracy between Bo and 77 
noted earlier, there is a clear minimum in and we therefore 
adopted the mean values of Bq and weighted by 1 jy^ as our 
best estimates of the magnetic-field parameters. 

We then made 35 simulations with the best-fitting values of 
Bo and 77 for each outer scale and evaluated the weighted ;t'^'s for 
the resulting cr^u profiles. All three values of A^ax we investi- 
gated give reasonable fits to the observed ctrm profile along the 
whole radio source. The fit for A^ax = 65 kpc is marginally bet- 
ter than for the other two values (x^^^ - 1 -8), consistent with the 
results of Sect. 17.31 below. In this case the central magnetic field 
strength is 3.5+1.2/iG and the break radius is 16+11 kpc. For 
the power spectrum with A^ax = 65 kpc and these best-fitting 
parameters, we also produced three-dimensional simulations at 
a resolution of 1.25 arcsec. Even though the fitting procedure is 
based only on the low-resolution data, this model also repro- 
duces the 1.25-arcsec profile very well (^^£(5=0. 7). Combining 
the values of for the two resolutions, using the 1.25-arcsec 
profile close to the core and the 5.5-profile at larger distances, 
we find ^2^^=1.8. 

Fig. [12] shows a comparison of the observed radial profiles 
for rms Faraday ctrm and (RM) with the synthetic ones derived 
for this model. The synthetic o-rm profile plotted in Fig.[T2]is the 
mean over 35 simulations, and may be compared directly with 
the observations. In contrast, the (RM) profile is derived from a 
single example realization. It is important to emphasize that the 
latter is one example of a random process, and is not expected 
to fit the observations; rather, we aim to compare the fluctuation 
amplitude as a function of position. 

The values of Bo and x^^^i for all of the three-dimensional 
simulations, together with 77 and for the single and double 
power-law profiles, respectively, are summarized in Table|6] 



7.3. The outer scale oi magnetic-field fluctuations 

The theoretical RM structure function for uniform field strength, 
density and path length and a power spectrum with a low- 
frequency cut-off should asymptotically approach a constant 
value (2cr^^ for a large enough averaging region) at separations 
~ Amax- The observed RM structure function of the whole source 
is heavily modified from this theoretical one by the scaling of 
the electron gas density and magnetic field at large separations, 
which acts to suppress power on large spatial scales. In Sect.|6] 
we therefore limited the study of the structure function to sub- 
regions of 3C 449 in which uniformity of field strength, density 
and path length (and therefore of the power-spectrum amplitude) 
is a reasonable assumption, inevitably limiting our ability to con- 
strain the power spectrum on the largest scales. 

Now that we have an adequate model for the variation of 
ne{r){B^(r)y^^ with radius (Eq.fTTTl. we can correct for it to de- 
rive what we call the pseudo-structure function - that is the struc- 
ture function for a power-spectrum amplitude which is constant 
over the source. This can be compared directly with the structure 
functions derived from the Hankel transform of the power spec- 



14 



D. Guidetti et al.: Structure of the magnetoionic medium around the FR Class I radio galaxy 3C449 



5.5— arcsec RM 


1 1 1 1 1 1 1 1 1 

A^g^ = 205 kpc 


- 






o 




o 






Bo=2.8 + /-0.5 fj.G 


(a) 


7J = 0.1+/-0.1 



1.25-arcsec RM 



A^^^=S05 kpc 



(b) 



B„=4.1+/-l.l /iG 
77 = 0.8 + /-0.4 



Fig. 9. (a) and (b): Distributions of the best-fitting values of Bq and rj from 35 sets of simulations, each covering ranges of 0.5 - 
lOyuG in Bo and - 2 in 77, at 5.5 and 1.25 arcsec respectively. The sizes of the circles are proportional to for the fit and the blue 
crosses represent the means of the distributions weighted by 1 /x^- The plot shows the expected degeneracy between Bq and 77. 
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Fig. 10. (a): Observed and synthetic radial profiles for rms Faraday ctrm at 5.5 arcsec as functions of the projected distance from the 
radio source centre. The outer scale is Amax=205 kpc. The black points represent the data with vertical bars corresponding to the 
rms error of the RM fit. The magenta triangles represent the mean values from 35 simulated profiles and the vertical bars are the 
rms scatter in these profiles due to sampling, (b) as (a) but at 1.25 arcsec. 



Table 6. Results from the three-dimensional fits at both the angular resolutions of 1.25 and 5.5 arcsec. 



FWHM 






single r] 






broken 77 


(arcsec) 


(kpc) 


So (mG) 


1 




So (tJG) 


r,„(kpc) xia 


1.25 


205 


4.1±1.1 


0.8±0.4 


0.48 






5.50 


205 


2.8±0.5 


0.1±0.1 


2.2 


3.5+0.7 


17±9 1.9 


5.50 


65 








3.5+1.2 


16±11 1.8 


5.50 


20 








3.5±0.8 


11+8 2.1 



trum. To evaluate the pseudo-structure function, we divided the 
observed 5.5-arcsec RM image by the function: 



neirYBirfdl 



1/2 



(12) 



(EnBlin & Vogt 2003), where the radial variations of and B are 
those of the best-fitting model (Sect.|7| and the upper integration 
limit L has a length of 10 times the core radius. The integral 
was normalized to unity at the position of the radio core: this is 
equivalent to fixing the field and gas density at their maximum 
values and holding them constant everywhere in the group. The 
normalization of the pseudo-structure function should then be 
quite close to that for the two central jet regions. 



The pseudo-structure function is shown in Fig.[T4f a) together 
with the predictions for the BPL power spectra with Amax=205, 
65 and 20 kpc. As expected, the normalization of the pseudo- 
structure function is consistent with that of the jets (Fig. [8^ and 
b). 

The comparison between the synthetic and observed pseudo- 
structure functions indicates firstly that they agree very well at 
small separations, independent of the value of Amax- This con- 
firms that the best BPL power spectrum found from a combined 
fit to all six sub-regions is a very good fit over the entire source. 
Secondly, despite the poor sampling on very large scales, the 
asymptotic values of the predicted structure functions for the 
three values of A^ax are sufficiently different from each other 
that we can determine an approximate outer scale. The model 
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Fig. 11. Distribution of the best-fitting values of Bo and from 35 sets of simulations, each covering ranges of 0.5 - lO/zG in 
Bo and - 50 kpc in r^. The sizes of the circles are proportional to and the blue cross represents the means of the distribution 
weighted by llx^- The plot shows the degeneracy between Bq and described in the text. 
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Fig. 12. Comparison between observed and synthetic profiles of rms and mean Faraday rotation at resolutions of 5.5 arcsec (a and 
b) and 1.25 arcsec (c and d). The synthetic profiles are derived from the best-fitting model with A^ax = 65 kpc. The black points 
represent the data with vertical bars corresponding to the rms fitting error, (a) and (c): Profiles of ctrm . The magenta triangles 
represent the mean values from 35 simulations at 5.5-arcsec resolution and the vertical bars are the rms scatter in these profiles due 
to sampling, (b) and (d): Profiles of (RM) derived from single example realizations at 5.5 and 1.25-arcsec resolution. 



with Aniax=65 kpc gives the best representation of the data. The 
fit is within the estimated errors except for a marginal discrep- 
ancy at the very largest (and therefore poorly sampled) separa- 
tions. That with Ainax=20kpc is inconsistent with the observed 
pseudo-structure function for any separation ^20 arcsec (^6 kpc) 
where the sampling is still very good, and is firmly excluded. 
The model with An,ax=205 kpc has slightly, but significantly too 
much power on large scales. 



We emphasize that our estimate of the outer scale of the RM 
fluctuations is essentially independent of the functional form as- 
sumed for the variation of the field strength with radius in the 
central region, which affects the structure function only for small 
separations. Our results are almost identical if we fit the field- 
strength variation with either the profile of Eq.[TO] (with ?; ^ 0) 
or that of Eq.fTTI 

As for the structure functions of individual regions, the 
pseudo-structure function at large separations is clearly affected 




Fig. 13. Comparison of observed and representative synthetic distribution of Faraday RM at 1 .25 arcsec and 5.5 arcsec. The synthetic 
images have been produced for the best-fitting model with Aniax=65 kpc. The colour scale is the same for all displays. 
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by poor sampling: this increases the errors but does not produce 
any bias in the derived values. At large radii, however, the in- 
tegral in equation [12] becomes small, so the noise on the RM 
image is amplified. This is a potential source of error, and we 
have therefore checked our results using numerical simulations. 
We calculated the mean and rms structure functions for sets of 
realizations of RM images generated using the FARADAY code, 
as in Sect. 17. II for different values of A^ax- These structure func- 
tions are plotted in Figs. [T4T b) and (c). 

The main difference between the model structure functions 
derived from simulations and the pseudo-structure functions de- 
scribed earlier is that the former show a steep decline in power on 
large scales in place of a plateau. This occurs because the smooth 
fall-off in density and magnetic field strength with distance from 
the nucleus suppresses the fluctuations in RM on large scales. 

The results of the simulations confirm our analysis using the 
pseudo-structure function. The mean model structure function 
with Amax = 65 kpc again fits the data very well, except for a 
marginal discrepancy at the largest scales. In view of the devi- 
ations from spherical symmetry on large scales evident in the 
X-ray emission surrounding 3C449 (Fig.[D^ we do not regard 
this as a significant effect. 

In order to check that the best-fitting density and field 
model also reproduces the data at small separations, we repeated 
the analysis at 1.25-arcsec resolution. The observed pseudo- 
structure function is shown in Fig. 14(d), together with the the 
predictions for the BPL power spectrum with Ainax= 205, 65 and 
20 kpc. The observed structure function is compared with the 
mean from 35 simulations with Amax= 65 kpc in Fig. 14(e). In 
both cases, the agreement for An,ax= 65 kpc is excellent. 

In Fig. [131 example realizations of this model with the best- 
fitting field variation are shown for resolutions of 1.25 and 
5.5 arcsec alongside the observed RM images. 

8. Summary and comparison with otiier sources 

8.1. Summary 

In this work we have studied the structure of the magnetic field 
associated with the ionized medium around the radio galaxy 
3C449. We have analysed images of linearly polarized emission 
with resolutions of 1.25 arcsec and 5.5 arcsec FHWM at seven 
frequencies between 1.365 and 8.385 GHz, and produced im- 
ages of degree of polarization and rotation measure. The RM 
images at both the angular resolutions show patchy and random 
structures. In order to study the spatial statistics of the magnetic 
field, we used a structure-function analysis and performed two- 
and three-dimensional RM simulations. We can summarize the 
results as follows. 

1 . The absence of deviations from rotation over a wide range 
of polarization position angle implies that a pure foreground 
Faraday screen with no mixing of radio-emitting and thermal 
electrons is a good approximation for 3C449 (Sect.|4|. 

2. The dependence of the degree of polarization on wavelength 
is well fitted by a Burn law. This is also consistent with 
pure foreground rotation, with the residual depolarization 
observed at the higher resolution being due to unresolved 
RM fluctuations across the beam (Sect. [5]). There is no evi- 
dence for detailed coiTelation of radio-source structure with 
either RM or depolarization. 

3. There is no obvious anisotropy in the RM distribution, con- 
sistent with our assumption that the magnetic field is an 
isotropic, Gaussian random variable. 



4. Our best estimate for the Galactic contribution to the RM of 
3C449 is a constant value of -160.7radm"^ (Sect.l^Ib. 

5. The RM structure functions for six different regions of the 
source are consistent with the hypothesis that only the am- 
plitude of the RM power spectrum varies across the source. 

6. A broken power-law spectrum of the form given in Eq.|9] 
with q\= 2.07, ^h=2.98, /b=0.031 arcsec"' and /^ax = 
1 .67 arcsec"' (corresponding to a spatial scale Amin=0.2 kpc) 
is consistent with the observed structure functions and depo- 
larizations for all six regions. No single power law provides 
a good fit to all of the structure functions. 

7. The high-frequency cut-off in the power spectrum is required 
to model the depolarization data. 

8. The profiles of ctrm strongly suggest that most of the fluc- 
tuating component of RM is associated with the intra-group 
gas, whose core radius is comparable with the characteristic 
scale of the profile (Sect. 14.11 ). The symmetry of the profile 
is consistent with the idea that the radio source axis is close 
to the plane of the sky. 

9. We therefore simulated the RM distributions expected for an 
isotropic, random magnetic field in the hot plasma surround- 
ing 3C 449, assuming the density model derived by Croston 
et al. (2008). 

10. These three-dimensional simulations show that the depen- 
dence of magnetic field on density is best modelled by a 
broken power-law function with B{r) oc rieir) close to the 
nucleus and B{r) ^ constant at larger distances. 

1 1 . With this density model, our best estimate of the central mag- 
netic field strength is Bo - 3.5 + 1.2yuG. 

12. Assuming these variations of density and field strength with 
radius, a structure-function analysis can be used to esti- 
mate the outer scale A^ax of the magnetic-field fluctuations. 
We find excellent agreement for A^ax ~ 65 kpc (fmin = 
0.0053 arcsec"'). 

8.2. Comparison with otiier sources 

Our results are qualitatively similar to those of Laing et al. 
(2008) on 3C 3 1 . The maximum RM fluctuation amplitudes 
are similar in the two sources, as are their environments. For 
spherically-symmetric gas density models, the central magnetic 
fields are almost the same: Bq ^ 2.8yuG for 3C31 and 3.5/iG 
for 3C449. Both results are consistent with the idea that the 
RM fluctuation amplitude in galaxy groups and clusters scales 
roughly linearly with density, ranging from a few rad m"^ in the 
much sparsest environments (e.g. NGC 315; Laing et al. 2006), 
through intermediate values 30 - 100 rad m"^ in rich groups 
such as 3C 3 1 and 3C 449 to ~ 10"^ rad m"^ in the centres of clus- 
ters with cool cores. 

The RM distribution of 3C 3 1 is asymmetrical, the northern 
(approaching) side of the source showing a much lower fluctu- 
ation amplitude, consistent with the inclination of 50° esti- 
mated by Laing & Bridle (2002). Detailed modelling of the RM 
profile led Laing et al. (2008) to suggest that there is a cavity 
in the X-ray gas, but this would have to be significantly larger 
than the observed extent of the radio lobes and is, as yet, un- 
detected in X-ray observations. A broken power-law scaling of 
magnetic field with density, similar to that found for 3C 449 in 
the present paper, would also improve the fit to the ctrivi profile 
for 3C 3 1 ; alternatively, the effects of cavities around the inner 
lobes and spurs of 3C449 might be significant. Deeper X-ray 
observations of both sources are needed to resolve this issue. 
In neither case is the magnetic field dynamically important: for 
3C 449 we find that the ratio of the thermal and magnetic-field 
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Fig. 14. (a): Comparison of the observed pseudo-structure function of the 5.5-arcsec RM image as described in the text (points) 
with the predictions of the BPL model (curves, derived using a Hankel transform). The predicted curves are for Amax=205 (blue 
dotted), 65 (continuous magenta) and 20 kpc (green dashed), (b) and (c): structure functions of the observed (points) and synthetic 
5.5-arcsec RM images produced with the An,ax=205,20kpc and An,ax=65kpc, respectively, (d) and (e) as (a) and (c) respectively, 
but at 1 .25-arcsec resolution. The error bars are the rms from 35 structure functions at a given A^ax, in (a) and (d) are representative 
error bars from the model with Amax=65 kpc. 



pressures is ^30 at the nucleus and ^400 at the core radius of 
the group gas, = 19 kpc. The magnetic field is therefore not 
dynamically important, as in 3C 3 1 . 

The magnetic-field power spectrum in both sources can be 
fit by a broken power-law form. The low-frequency slopes are 
2. 1 and 2.3 for 3C 449 and 3C 3 1 respectively. In both cases, the 
power spectrum steepens at higher spatial frequencies, but for 
3C31 a Kolmogorov index (11/3) provides a good fit, whereas 
we find that the depolarization data for 3C449 require a cut- 
off below a scale of 0.2 kpc and a high-frequency slope of 3.0. 
The break scales are »5 kpc for 3C 31 and a;l 1 kpc for 3C449. 
It is important to note that the simple parametrized form of the 
power spectrum is not unique, and that a smoothly curved func- 
tion would fit the data equally well. 

The gas-density structure on large scales in the 3C 3 1 group 
is uncertain, so Laing et al. (2008) could only give a rough 
lower limit to the outer scale of magnetic-field fluctuations, 
Amax ~ 70 kpc. For 3C449, we find A^ax ~ 65 kpc. The pro- 
jected distance between 3C449 and its nearest neighbour is 
~33kpc (Birkinshaw et al. 1981), similar to the scale on which 
the jets first bend through large angles (Fig.[T]i. As in 3C 31, it is 
plausible that the outer scale of magnetic-field fluctuations is set 
by interactions with companion galaxies in the group. 
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